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Abstract 



We tabulate bounds on the optimal number of mutually unbiased bases in M . For most dimen- 
sions d, it can be shown with relatively simple methods that either there are no real orthonormal 
bases that are mutually unbiased or the optimal number is at most either 2 or 3. We discuss the 
limitations of these methods when applied to all dimensions, shedding some light on the difficulty 
of obtaining tight bounds for the remaining dimensions that have the form d = 16n^, where n can 
be any number. We additionally give a simpler, alternative proof that there can be at most d/2 + 1 
real mutually unbiased bases in dimension d instead of invoking the known results on extremal 
Euchdean hue sets by Cameron and Seidel, Delsarte, and Calderbank et al. 

Keywords: Quantum Information Processing, Quantum Computing, Hadamard Matrices, Euclidean 
Line Sets, Hadamard Conjecture. 

1 Introduction 

Two orthonormal bases B and B' of the Hilbert space are called mutually unbiased if and only if 
the modulus of the standard complex inner product is 



for all \(f)) £ B and all \ip) £ B' . The problem of determining bounds on the maximum number M^d of 
bases over that are mutually unbiased is an important open problem |15| which has received much 
attention |E1 El 121 EI • We refer the reader to e.g. [T] for an overview of known bounds. 

The most common application of MUBs is found in quantum cryptography where MUBs are the 
quantum states used in most QKD protocols [3 01 El 15 ■ 

The MUB problem has remained open for over 20 years. For example, it is not even known whether: 
• for all M there is a d' such that M^d > M for all d> d'l In this manuscript, we consider the problem 
of determining M^d the maximum number of mutually unbiased bases over R'^ {real MUBs), defined 
analogous to the above. In the real case, the problem appears to have a somewhat different character: 
elementary arguments give tight bounds on almost all d] and the answer to the question • is easily seen 
to be negative. Additionally, we are able to directly construct real MUBs which cannot be extended 
to an optimal set even using complex bases; which means in general, maximal sets of MUBs are not 
necessarily optimal. 

* corresponding author wocjan@cs.caltech.edu 
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Table 1: Here s is any positive odd integer, and i is any non-negative integer; MOLS{r) denotes the 
maximum number of mutually orthogonal Latin squares of order r. Recall that a Hadamard matrix 
(abbreviated by HM) is a square matrix of size d with entries in {—1, 1} with pairwise orthogonal 
columns. 



Obtaining real MUBs is a special case of the geometric problem of obtaining "extremally or uni- 
formly distributed" collections of Euclidean line sets JHI with arbitrary, prescribed angles. Observe 
that a collection of MUBs is exactly a collection of orthogonal sets of lines passing through the origin 
such that any one of the following equivalent statements holds: (i) the angle between any pair of 
lines taken from different sets is the same; (ii) the minimum angle over any pair of lines taken from 
different sets is maximized; (iii) the maximum angle over any pair of lines taken from different sets 
is minimized. (This formulation of MUBs applies to any Hilbert space when "angle" is replaced by 
"magnitude of inner product.") 

In this paper, we show upper and lower bounds on the number of real MUBs in various dimensions 
using elementary methods or known results |17l llUj . These results are summarized in Tabled The 
paper is organized as follows: Section [21 introduces a simple proof technique which we use throughout 
this paper. We use a simple technique to find the number of real MUBs in dimensions which are not 
divisible by four {M^d = 1) and non-square dimensions {M^d < 2). The remaining dimensions can be 
written as d = 4*s^ with s odd and z > 1. We consider such square dimensions 4*s^ in Section |31 where 
we find an upper bound (Mjgd < 3) in the case of i = 1, and we use a previous result [1,7^ to show a 
lower bound in the case of i > 1. Section 0] gives a simpler, alternate proof for a previously known 
upper bound {M^d < d/2 + 1). 

Finally, in Section El we consider the question of whether all sets of real MUBs can be extended 
into an optimal set of (complex or real) MUBs and find that the answer is negative. 



2 Most (i's admit at most 1 or 2 real MUBs 

In this section, we consider non-square dimensions and dimensions not divisible by four. The proofs 
use an elementary technique which we will make use of in later sections. We will repeatedly use the 
following standard canonical form for real MUBs. 

Observation 2.1. If Bi, B2, . . . , B fji are a collection of real MUBs, without loss of generality, we can 
always choose Bi to be the standard basis (simply multiply all matrices by B^^ from the left). Putting 
the elements of Bi as columns, we get the identity matrix 1. By Equation (1), all the other Bi are 
orthogonal and unbiased to Bi if and only if their basis elements are the columns of a Hadamard 
matrix scaled by and Bi,Bj are pairwise unbiased if and only if B^^Bj = BjBj is a Hadamard 

matrix scaled by 
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Proposition 2.2. If A )(d, then M^a = 1. 

Proof. Let Bi, B2 be two real MUBs in M"^. Using the canonical form of Observation 12.11 for Bi and 
B2, it is a simple folklore fact that "d must therefore be 2 or a multiple of 4 since otherwise Hadamard 
matrix of size d cannot exist." (The famous Hadamard conjecture ^2] is that for every multiple of 4, 
there exists in fact a Hadamard). 

We prove this simple folklore fact here, since many of the arguments presented in this note are 
elementary variants of this proof. This proof is similar to the one found in |16l page 74]. We show the 
stronger statement that no more than 3 orthogonal vectors in { — 1, l}'^ can exist unless d is a multiple 
of 4. 

Assume there are three vectors vi,V2,V3 E {—1,1}'^ that are orthogonal to each other. Observe 
that permuting the entries of the three vectors according to any permutation does not change the 
inner product of any pair of them. Also note that changing the sign of the entries by multiplying the 
three vectors entry-wise with any column vector does not change the inner product of any pair, nor 
does it change the magnitude of any entry. 

We use these two facts to assume without loss of generality, that the first vector vi is the all ones 
vector (+,+, . . . ,+). The second vector V2 should be orthogonal to the first, so again without loss 
of generality it can be chosen as the vector (+, +,...,+,—,—,..., — ), with the first half +'s and the 
second half — 's. This requires d to be even and completes the Hadamard conjecture for d = 2. The 
third vector V3 should be orthogonal to the first and the second, hence without loss of generality, it 
can be chosen to be of the form (+, ...,+,—,...,—,+,...,+,—,..., — ), with 4 blocks of +'s and — 's 
with alternating signs. This is due to the ability to permute the entries in any way we like. Any 
permutation that fixes vi and V2 is allowed, so we permute such that there are ?n.++ values of + 

followed by values of — in the first half of the vector. In the second half of the vector there are 

m |_ values of + followed by m values of — . This gives rise to 4 independent equations: 

d 

d 

m L + m = — 

+ 2 

{vi , V3) = m++ — + m ^. — m = 

{^2,^3) = — — m ^ + m = 

The above equations give rise to the unique solution m±± = d/A and consequently A\d. □ 

This shows that there are two real orthonormal bases are that are mutually unbiased whenever 
4 / d. Next we consider non-square dimensions of the form d = An. 

Proposition 2.3. If A\d, hut d is not a square, then Mjgd < 2. In this case, M^d = 2 if and only if a 
Hadamard of order d exists. 

Proof. Assume the contrary and let Bi, B2 and B3 three such orthonormal bases that are mutually 
unbiased. We also use -61,-62 and B3 to denote the corresponding real matrices in the standard 
canonical form of Observation 12.11 

As in the proof of the simple folklore fact and Proposition 12.21 given above, we can always achieve 
without loss of generality that the first column of B2 is -^(1, 1, ... , 1). This is done as follows. If the 

ith entry of B2 is ^ then multiply the ith. row of B2 and the ith row of -B3 by — 1. Clearly, this does 
not change the absolute values of the inner products between the columns vectors of Bi, B2 and B^, 
nor does it change the orthogonality of any of them. 
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Now, let -^(si, S2, ■ ■ ■ , Sd) (with Si = ±1) be the first column of B3. Eq. (1) applied to the first 
column of B2 and the first column of B3 implies that 

d 

\Y,Si\ = Vd. (2) 

1=1 

The absolute value of the sum is clearly a natural number, whereas \/d is irrational if d is not a square. 
Therefore, there cannot exist a third orthonormal real basis B3 that is mutually unbiased to Bi and 
B2. 

This shows that the maximum number of real MUBs is at most 2 for dimensions d that are not 
squares. However, if there is a Hadamard matrix of order d, then in fact, B2 exists, and we have 
exactly 2 real MUBs. □ 

In summary, we have handled the case of all dimensions not divisible by 4 by showing that there 
are not even two MUBs. In non-square dimensions which are divisible by four, there are at most two 
MUBs. In the next section, we will consider all other dimensions. 



3 Square dimensions divisible by 4 

We have so far shown that unless d = Am? for some positive integer m, there are at most 2 real MUBs, 
such dimensions can be written as d = 4*5^ with s odd. First we consider the case of z = 1 and then 
consider the case of i > 1. 

Proposition 3.1. If d = 45^ for an odd positive integer s, the M^d < 3; furthermore M^d > 2 
provided that a Hadamard matrix of order d exists. 

Proof. We show a stronger Lemma 13.31 that involves mutually unbiased lattice lines (as opposed to 
entire bases) . Using this lemma and the canonical form of Observation 12.11 it follows that only 2 
further MUBs exist that are unbiased to the standard basis. □ 

Definition 3.2. A lattice line in is a line passing through the origin and some point v in {—1, l}*^. 
We denote the line by either of the vectors v or —v. Two lattice lines v and w are mutually unbiased 
if \{v,w)\ = where (•,•) is the standard inner product in . We denote by L^d the maximum 

number of mutually unbiased lattice lines that can be found in M.'^. 
Lemma 3.3. If d = 4.^ for an odd positive integer s, then L^d < 2. 

Proof. Assume to the contrary that there are 3 such lines represented by the vectors vi,V2,Vd S 
{ — 1,1}'^. Consider the vectors corresponding to those lines. As in the proof of the simple folklore 
fact and Proposition 12.21 we assume that vi = (+,+,...,-1-) is the all-one-vector and we consider 

a partition of these vectors into 4 blocks with lengths m++ , , m |- , and m We obtain the 

following equations: 

{vi,Vi) = m++ -|- m-i -|- m ^ -|- m__ = 4.^ 

(wi, 173) = — -|- m ^ — m__ = ib2s 

{vi,V2) = + m_| — m ^ — m__ = ib2s 

{'^2,^3) = — — m ^ -|- m__ = ib2s 
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One can solve for m++ by adding all the equations together to get: 



m_|__ 



2 s(±l±l±l) 
s + — 7: 



Since the value in the parenthesis can only add up to an odd number (—3, —1, 1, 3), and due to the 
fact that s is odd and the product of two odds is odd, '^(^^^^^^^ cannot be an integer. However, if 
m++ is not an integer, we have a contradiction, proving the result. □ 

Hence we see that if d = 4*s^ with s odd and i = 1, there are at most 3 real MUBs. Now we 
consider the other half of the remaining cases, i.e., where i > 1. In this case, an earlier construction 
of MUBs using Latin squares immediately gives a lower bound. 

Proposition 3.4. |.?7| / If d = 4*s^, where s is any positive integer, then M^d > MOLS{T- s) + 2, 
provided that there exists a Hadamard matrix of order 2*s, where MOLS{m) denotes the maximum 
number of mutually orthogonal Latin squares of order m. 

This construction is contained in Appendix Since there is a Hadamard of order 2, the above 
construction also works in dimension 4 and is optimal in that case. We obtain the following three 
mutually unbiased bases: 
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So far, we have seen that the upper bound on the number of real MUBs in dimension d is constant 
(either 1,2 or 3) unless d = 4*s'^ with s odd and i > 1. In the latter case, we found a lower bound 
on the number of real MUBs following from the construction based on mutually orthogonal Latin 
squares. In the next section we give an upper bound on the number of real MUBs and see that the 
upper bound is tight in dimension d = 4* . 



4 A general upper bound which is tight for (i = 4* 

In this section we consider general dimensions and give a new proof that M^d < d/2 + 1. Then we 
cite a construction for d = 4* which obtains the upper bound, and thus we see that the bound is tight 
in the case of d = 4*. 

Proposition 4.1. [Wjl We have the upper bound M^d < d/2 + 1 for any d. 
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Proposition 14. II is a special case of the type of bound proved by [S] (using Jacobi polynomials) and 
jlUj (using tensor algebra first principles), where they consider the maximum number of lines that can 
be packed into such that the angle between any pair of lines is one of two specified angles. 

Next, we give an alternative, simpler proof of this upper bound. This is based on an adaptation 
to M!^ - of a result of |2j which equates the MUBs over to so-called commuting classes of unitary 
matrices of size d x d. 

Theorem 4.2. ^2] There exists a set of m MUBs in if and only if there are m classes Ci, . . . , Cm 
of the following properties: 

• each class Cj consists of d commuting matrices, 

• any two classes have only the identity matrix in common, and 

• all matrices in Ci, . . . , Cm are pairwise orthogonal with respect to the trace inner product. 

The corresponding MUBs Bi,. . . ,Bm are the common eigenvectors of the matrices of the com- 
muting classes Ci, . . . , Cm, respectively. 

Corollary 4.3. There exists a set of m MUBs in M'^ if and only if there are m classes Ci, . . . ,Cm 
with properties as in Theorem \4-^ with the additional property that all matrices are real symmetric. 

Proof. Since the eigenvectors of real symmetric matrices are real, it follows from Theorem 14.21 that 
there are m real MUBs under the assumptions of the corollary. This proves one direction of the 
equivalence. 

For the other direction, we show how to construct m commuting classes Ci, . . . , Cm of real sym- 
metric matrices satisfying the properties of Theorem 14.21 starting from m MUBs. 
Let the j'th MUB be given by 

B, = {\4),...M)}- 

Assuming that m > 2, we can construct at least one real Hadamard 

H := VdB\B2 , 

where we use Bi and B2 to also denote the unitary matrices whose column vectors are given by basis 
vectors of Bi and B2, respectively. We may assume without loss of generality that the Bi is the 
standard basis and that the first vector of B2 is given by the normalized. Then the first row of H is 
the all-one- vector. Denote the entries of H by htk- 

Then we define the jth commuting classes Cj := Uj^2, ■ ■ ■ , f^j.d}) where the matrices within 

the class are given by 

d 

k=l 

for t = 1, . . . ,d and j = 1, . . . , m. 

Observe that the matrix Uj^t is diagonal with respect to the basis Bj and that its eigenvalues are 
given by the entries of the tth row of the real Hadamard matrix H. Therefore, theses matrices are 
not only unitary but also real symmetric. Clearly the matrices within each commuting class commute 
because they are diagonal with respect to the same basis. The only matrix that is contained in all 
commuting classes is / = Uj^i. It remains to show that they are all orthogonal with respect to the 
trace inner product. 

Observe that all matrices within a commuting class are orthogonal because the rows of H are 
orthogonal. This implies that all matrices not equal to the identity matrix are traceless. 
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We can compute the inner product directly: 



tr{ulUj,s) = tr{Ui,rUj,s) 



x=l y=l 
^ d d 

x=l y=l 

1 

d 



because one of the matrices is not the identity matrix. This completes the proof. □ 

Now we are ready to give the alternative proof of Proposition 14.11 

Proof, (of Proposition 14. If there are m real MUBs, then according to Corollary 14.31 we can con- 
struct m classes of commuting real symmetric unitaries. Each class contains d — 1 traceless matrices. 
Additionally, we can consider the identity matrix, which is also a real symmetric unitary. All of these 
m{d — 1) + 1 matrices are orthogonal. The space of real symmetric matrices is a vector space of 
dimension d{d + l)/2. The span of our m{d — 1) + 1 matrices is a subspace of all real symmetric 
matrices, thus: 



m{d + 1 < 



d{d+l) 



2 

d 

This constitutes a new proof for the upper bound of real MUBs in dimension d. □ 

A construction of d/2 + 1 line sets that correspond to real MUBs is given in for the special 

case of d being a power of 4. The general upper bound of Proposition 14. II shows that construction is 
optimal. 

In the next section we consider the question of extending sets of MUBs to try to reach the upper 
bounds of this section. Unfortunately, we will see that such extensions in general are not optimal. 



5 Greedy methods do not work 

In this section, we consider the efficacy of greedy methods of constructing real MUBs. We ask the 
question: Starting with an arbitrary set of MUBs or mutually unbiased lattice lines, can one extend 
the set to End an optimal number of MUBs or mutually unbiased lattice lines? Put another way: can 
we take any MUB that comes our way, or must we be careful to choose sets that fit well together. 

We consider the question: can the method of Proposition 13.11 and Lemma 13.31 be extended to 
show an upper bound on the number of real MUBs for d = 4*s^ (to complement the lower bound of 
Proposition 13. 4|) ? 

We could continue, without loss of generality, to partition j vectors into 2^ blocks and write down 
equations that they must satisfy in order to represent mutually unbiased lattice lines. Showing the 
conditions on j and d under which such systems do not have solutions would then give us an upper 
bound on L^d (and M^d) for specific dimensions d. However, in this approach, we end up with a 
family of linear Diophantine systems one for each j and d, with inequalities and equalities, since the 
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solutions rriij (to the sizes of the blocks in the partition) should be non-negative integers. Showing 
the conditions on j and d for which such a system does not have a solution does not appear to be 
tractable. 

Our approach is instead the following. We construct a particular set of 2* lattice lines in dimension 
d = 4ts^ using the 2* x 2* Sylvester Hadamard (I2j. We then show that this particular set cannot 
be extended by even one lattice line. If any set can be extended into an optimal set, then Ljgd = 2* 
for d = 4*s^. However, this result is not true, because we later give a construction that gives L^d > 
2%2\s^ - s) for some cases of i, which shows that not every set can be extended into an optimal set. 

Let H be the Sylvester Hadamard matrix of order 2* jl2j and denote its entries by Ski for /c, I = 
0,...,N — 1, where N = 2'. Let 'Wo,wi,...,W]\f-i be all-one-vectors of lengths 5o, 6i, . . . , ^at-i, 
respectively. The lengths will be determined later. The vectors vo,vi, . . . ,vn-i have the following 
block structure: 

Vk ■= {Sk,0 Wo I Sk,lWi I ... I Sk,N-l WN-lY , 

where | means that we concatenate the sub- vectors Sij Wj and t denotes the transposition. The condi- 
tions for the vectors to for a collections of lattice lines are given by the following equations 



{Vk,Vl) 
{Vk,Vl) 



2's 



for < i < j < iV — 1. It may appear that there are 2* -|- 2*(2* — l)/2 equations. But that is not so. 
Many of the equations are equivalent to others. This is because the columns of the Sylvester Hadamard 
matrix form a group^ under point-wise multiplication. So there are only 2* equations, rather than the 
much larger number that appears above. More precisely, if © denotes the group operation, then we 
have 

{Vk,Vl) = {vo,Vk(Bl) ■ 

So, in fact we see the only equations we have to consider are of the form 

N~l 



{vo,vo) = Z^bj = 4*s , {vo,Vk) 

j=0 

for 1<A;<A^— 1. In matrix notation, this means 



j=0 



±2's. 



( 



H 



bo 
bi 



\ 



( 4*s2 
±2*s 



By multiplying both sides by ff*/2* we obtain the solution 



bo \ 
bi 



\ bN-i J 



( 2's^ \ 
±s 

\±s ) 



^This group is isomorphic to x ■ ■ ■ x (the i-fold direct product of the group Z2). The multiphcation of the fcth 
and Ith column gives mth column, where m is obtained by XORING the binary numbers corresponding to k and I and 
converting the resulting binary number to decimal notation. We denote this operation by m = fc © L 
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and it is clear that the resulting lengths bo,bi, . . . , b^^i are integers. This completes the proof that 
there are at least 2* lattice lines in dimension d = 4*s^. 

Next, we show that such a collection of Sylvester lattice lines cannot be extended. Let vq, . . . , vjy-i 
be a collection of AT = 2* such lines in dimension d = 4*s^. Assume there is vector w that is mutually 
unbiased to the vectors vq,. . . ,vn-i. We may assume without loss of generality that in the bj block 



6^ + 6^. 



This 



of w consists of a sub-block of +l's of length b'j and a sub-block of — I's, where bj 
follows from the fact that permuting the entries of the vectors within the jth block according to the 
same permutation does not change the vectors vq, . . . ,f7v-i and does not change the inner product 
between w and vq, . . . , vj^-i. This leads to the following equation system: 



/ ^0 + &o \ 



(I2 H2 



h' -I- h" 



7V-1 



b'n 



N-1 



h" 
"0 




h" 



b'N-1 
\ b"N-i J 



= H2i+i (60 , 60 , . . . , b'j^_ I , b%_ 1 )* 
= P(4*s2,±2^s,...,±2*s)*, 

where P is a permutation matrix. Multiplying by i/gi+i/^*^^ on both sides we obtain 



,■ 1 n S(±l±l 



±1 



sk 



for some odd number k corresponding to the sum in the parenthesis. But since s is also odd the 
resulting number 6q is not an integer. This is a contradiction completing the proof that the collection 
of 2* Sylvester lines cannot be extended. 

Unfortunately, it turns out that for d = 4*s^, our greedy construction gave 2* unbiased lattice lines, 
but we can do a lot better for i > 1 as shown below. 

Specifically, in the (generalized) construction of Lemma given above, while the partitioning 
and choices of the first three vectors (needed for i = 1) are without loss of generality, the choice of the 
4th and further vectors (needed for i > 2) is specific and loses generality. Therefore, non-extendibility 
does not imply optimality, except for i = 1. 

Proposition 5.1. We have L^d > d — \fd 'provided that there is a Hadamard matrix of order d — \fd. 

Proof. Assume that there is a Hadamard matrix H of order d — \fd. We obtain d — \fd lattice lines 
in dimension d by appending to the row vectors of H the all-one- vector of length \fd. 

In other words, if the /ifc is the /cth row vector of -ff, then the k lattice vector has the form 

= {hk,i, hk^2, ■■■ , \,(rf_v^)! 1) 1, • • • ) 1) ) 
It is obvious that these vectors are lattice lines. □ 



The above proposition along with many examples gives counter-examples to the greedy conjecture. 
For instance consider d = 4?. If any maximal set were optimal, then L^d < 2^. But since there is a 
Hadamard of dimension 12 {d — Vd), we see that L^d > 12. 

So far, we have shown that greedy methods for constructing mutually unbiased lattice lines do 
not work. Next, we consider the efficacy of greedy methods of constructing general (complex) MUBs. 
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Specifically, we give a negative answer to the question: is every maximal, or non-extendible set of 
MUBs necessarily optimal? 

We consider the general, complex version of the MUBs constructed in ProDosition l3.4l '17' based on 
Latin squares, we call these the Latin MUBs: note that this construction yields at most y/d+ 1 MUBs 
for dimensions of the form d = p^^ , where p is a prime and e is an integer, whereas in prime power 
dimensions it is known that there are in fact d+\ MUBs j2] . Nets and the Latin MUB construction 
are reviewed in Appendix We obtain restrictions on possible extensions to latin MUBs in general 
and show that the set of latin MUBs in d = 4 is maximal, i.e., not extendible, thus answering the 
above question in the negative. 

2 

Observation 5.2. The s{s + 1) incidence vectors of a {s + 1, s)-net span . 

Proof. For any pG{l,...,s^}, there is exactly one incidence vector in each parallel class that has 1 as 
entry at position p. Sum all these vectors. This gives the vector t;p = (1, . . . , 1, (s + 1), 1, . . . , 1), where 
the value (s + 1) appears at position p. Now the s^ vectors Vp span W because the pth standard basis 
vector Cp is simply: (l/((s + l)s^)((s^ + s — l)vp — Vg). □ 

Lemma 5.3. In dimension d = s^, any new basis unbiased to a {s + 1) collection of Latin MUBs 
corresponds to a Hadamard matrix. (Note the standard basis is not unbiased to any Latin MUB). 

Proof. Assume there is a (s + 1, s)-net and let us consider the {s + 1) Latin MUBs that are constructed 
based on this net. We show that any orthonormal basis of that is mutually unbiased to the Latin 
MUBs defines a generalized Hadamard matrix of order d, that is, the absolute values of all entries is 
1/Vd. 

Let us consider the ith incidence vector of the block b. This incidence vector gives rise to s 
basis vectors of the 6th Latin MUB. Denote the support of m^j by Si. By construction the s basis 
vectors contain a Hadamard matrix H of order s as a sub-matrix if we look at the entries at positions 
contained in Si. 

Let y = {yi, . . . , y^) be any vector of the new ONB and denote by yi = {yi^ , . . . ,yi^) the sub- vector 
of length s of y defined by the support Si = {ii, . . . ,is} C {I, . . . ,d}. Since y is mutually unbiased to 
the s basis vectors obtained from m^j, we have 

-^Hiji = u (3) 

where u = [ui, . . . , Ug) is some vector in with \uj\'^ = ^ for j = 1, . . . , s. Since the matrix is 
unitary, it follows that = We have = s/d = 1/s and consequently 

= i'^bi, ■ ■ ■ , \yd\y) = ^ ■ 

The above equation is true for all incidence vectors m^j and gives a system of linear equations for the 
squares of the absolute values 

Since the support of each incidence vector has cardinality s, one solution to the above system is 
by choosing all entries of y to have the same magnitude, that is, \yj\'^ = Due to Observation 15.21 
the system has full rank and thus this solution is always a unique solution. This complete the proof 
that the entries of all basis vectors of the new ONB have the same magnitude and consequently that 
corresponds to a generalized Hadamard matrix. □ 

Proposition 5.4. The 3 standard Latin MUBs in dimension d = 4 are unextendible. 
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Proof. Assume there is a vector y that is mutually unbiased to all vectors of the three Latin MUBs. 
It follows from Lemma 15.31 that all its entries have the same modulus. Without loss of generality we 
may assume that y = (1, a, b, c)/2 after factoring out a common phase factor. 

Writing the equalities |1 + = |1 + 6p = |a + 6p = 2 that follow from eq. Olit is easily seen that 
there is no such a and b satisfying these conditions. This shows that we cannot even find a single 
vector that is unbiased to the Latin MUBs for d = 4. This completes the proof that the three Latin 
MUBs cannot be extended. □ 

In dimension d = 4 if one allows complex MUBs, 5 MUBs may be found. However, if we start 
with the three real Latin MUBs, we cannot even find a fourth MUB, complex or real, unbiased to 
the previous three. Hence it is clear that a set of MUBs which cannot be extended is not necessarily 
optimal. 

6 Conclusion 

In this paper we have found upper and lower bounds on the number of real MUBs in all dimensions. 
In most dimensions, the problem is solved. Assuming the Hadamard conjecture is true, the only 
interesting dimensions left are d = 4*s'^, where s is odd. In the case of d = 4s^, assuming the 
Hadamard conjecture is true, the only question is: does M^d = 3 hold always, or only for certain 
values of s. Additionally, what is M^d for d = A^s^ with i > 1 and s odd? Assuming the Hadamard 
conjecture is true, this paper shows a lower bound valid for all i > 1. It would be interesting to find 
an example that exceeds that lower bound. 

The work of M.S. and M.T. supported in part by NSF Grants EIA 02-18435 and CCF 04-04116. 
The work of P.W. supported in part by NSF Grant EIA 00-86038 and by the NSA under ARO contract 
number W911NF-05-1-0294 and by the NSF under contract number PHY-0456720. 
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A Nets and Latin MUBs 

We quote the definition of a net, which is used in the Latin MUB construction jJTj- Note that nets are 
equivalent to orthogonal arrays and traversal designs [16] . 

Definition A.l (Net). 

Let {mil, • • • ) mis) ^21, ■ ■ ■ , ^2s, ■ ■ ■ , mfci) • • • j ^ks} be a collection of ks incidence vectors of size d = 
that are partitioned into k blocks where each block contains s incidence vectors. The incidence 
vectors are denoted by m;,j, where b = 1, . . . ,k identifies the block and i = I, . . . , s the vector within a 
block. If the incidence vectors satisfy the following conditions we say that they form a {k,s)-net. 

1. The supports of all vectors within one block are disjoint, i.e., 

^li mfoj = (4) 

for all 1 <b < k and all 1 < i ^ j < s. 

2. The intersection of any incidence vectors from two different blocks contains exactly one element, 
i.e., 

m^i T^cj = 1 (5) 

for all 1 <b ^ c < s and all 1 < i, j < s. 
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Note that our definition of (k, s)-nets is in accordance with the usual definition of nets in design 
theory 8, page 172]. 

Let m E {0, 1}"' be an incidence vector of Hamming weight s and h G C an arbitrary column 
vector. Then we define the embedding of h into controlled by m, denoted by h | m, to be the 
following vector in 



h I m := h[r] |jV 



(6) 



r=l 



where h[r] is the rth entry of the vector h, {ji,j2, ■ ■ ■ ,js} the support of m with the ordering ji < 
j2 < . . . < js and \ jr) the jVth standard basis vector of C^. A less formal way to define this vector is: 
the first non-zero entry of m is replaced by the first entry of h, the second non-zero entry of m by the 
second entry of h, etc. 

This operation is best illustrated by a simple example: 

/ 1 \ 



1 

G{0,1}^ h:=| w |gC3, hTm:= G 
\ uj" 




V 1 J 

Theorem A. 2 (Construction of MUBs). 

Let {mil, . . . , mis, m2i, • • • , m2s, . . . , nifci, . . . , nifcj,} a {k,s)-net and H an arbitrary generalized 
Hadamard matrix of size s. Then the k sets for b = 1, . . . ,k 



m 



/ 


1 


\ 











UJ 







































V 




/ 



:= |-^(h, t rubi) I / = 1,..., 



s , i = 1, . . . , s 



(7) 



are k mutually orthogonal bases for the Hilbert space C^. 
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